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Abstract 



We study a weakly coupled 3 dimensional tp'^ type model consisting of N real 
5^ I scalars 99* coupled to an abelian gauge field A'^ and one extra scalar field p. We 

argue that, below some scale A = 0{T), this is the effective field theory of a 3+1 
dimensional abelian Higgs model at a high temperature T. The effective theory is 
sufficient to study the nature of the phase transition in four dimensions. By intro- 
ducing an auxiliary field x we eliminate the explicit ip"^ term; the new Lagrangian 
allows for a simple computation of the dominant corrections to the effective scalar 
potential V^ff in the large limit. We study three cases: a) Qg^/X < 0(1), b) 
6g^/X ~ 0{N) and c) 6g^ /X ~ 0{N'3) where g and A are the 4d gauge and scalar 
couplings, respectively. For case a) which is the most thoroughly studied we find 
that the leading 0{N) + 0(1) result for 14// admits only a second order phase tran- 
sition. For the other cases we find that b) the leading 0(1) result for 14// admits a 
first order phase transition whose strength is independent of and c) the leading 
0(A^3) result for Veff admits only a second order phase transition - the 0(1) cor- 
rections to this can be interpreted as indicating a first order phase transition whose 
strength diminishes as N increases. 



1. Introduction. 



There has been much recent interest in the nature of the electroweak phase transition. 

In general, for a spontaneously broken 4d 0^ type model, one can attempt to compute the 
high temperature effective scalar potential Veff. At very high temperature the symmetry 
is believed to be restored, and at some temperature T = T2 the origin (|0| = 0) goes from 
being a local minimum to being a local maximum. If at T > T2 the origin is also the global 
minimum the theory admits only a second order phase transition which proceeds by a 
roll-over after the temperature drops below T2. On the other hand, if at some Ti > T2 the 
scalar potential has another minimum, degenerate in energy with the one at the origin, 
the theory can admit a first order phase transition. 

The determination of a reliable effective potential has proven difficult because, in 
generic models, ordinary perturbation theory in h has infrared divergence problems just 
at values of 0^ and T where the phase transition occurs. In a 0^ model with real 
scalars, Dolan and Jackiw [1] showed over 18 years ago that the one-loop potential has 
problems. They summed up an infinite class of diagrams, the "superdaisies" to obtain a 
reliable result for the effective mass at the origin. 

Giving a reliable estimate of Veff in the case of the standard model is much more 
involved. The typical approach [2,3,4,5] for the corrections from the bosonic sector of the 
standard model has been to use not tree-level propogators but propogators corrected by 
one-loop vacuum polarization effects in a one-loop calculation of the effective potential. 
This has the effect of summing up an infinite class of diagrams, the so-called "ring" di- 
agrams, with the correct combinatorics except for the two-loop graphs [4,5,6] which are 
overcounted. When only the leading dependent corrections are used for the vacuum 
polarizations the overcounting is not serious, but when field dependent vacuum polariza- 
tions are used it can result in dangerous T^l^l terms in the effective potential [4,5] which 
results in the effective mass blowing up as one approaches the origin. 

Assuming no such dangerous terms appear, previous studies have concentrated on 
determining the coefficient of the |0|^ term in their estimates of Veff. For small 0^, the 
effective potential can be expanded as {d > 0) 

Veff ^c + acj)^ -b\(l)\^ + (!(!)'' + ■■■. (1) 

The coefficient a is positive when T > T2 and negative when T < T2. In a general 
3-1-1 dimensional scalar and gauge system the coefficients a, b and d can be found in 
the ring approximation to Veff as follows [2,3,4,5]. In the background field evaluation 
of Veff one expands the action about quantum scalar fields 0* and quantum gauge fields 
A'^. Each quantum degree of freedom has a mass depending on the background field 0^. 
At zero external momentum, the leading vacuum polarization effects give additional 
contributions to the mass of the quantum scalars and the longitudinal gauge field, but not 
to the transverse gauge fields. Then, when T is sufficiently bigger than T2, the result of 
the polarization effects is that the scalar and longitudinal gauge boson loops give only 0^ 
and (j)^ contributions to eq. (1). In this limit, the transverse gauge fields still contribute 
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a 101^ term. It is the existence of this term in the context of the standard model which 
led the authors of [2,4] to conclude that the electroweak model with sufficiently small 
self-scalar coupling admits a (weak) first order phase transition. 

The analysis we have just described is however incomplete. For example, how are 
these results modified if one keeps the field dependent contributions to the vacuum polar- 
izations? As already mentioned, the simple ring diagram procedure can lead to spurious 
contributions to K//, so obviously some care is needed. One would also hke to 

investigate how the momentum dependence of the vacuum polarizations modifies the |0^| 
result. Finally, wc would like to know in what way the ring (or "daisy") diagram sum can 
be interpreted as a consistent approximation scheme. 

To address these questions and problems, we study as a toy model a 3 dimensional (p^ 
type model consisting of N real scalars coupled to an abelian gauge field and an 
extra scalar p, which we argue is the low energy effective field theory for a 3+1 dimensional 
Abelian Higgs model at high temperature. The computations in the effective 3d theory are 
however much easier than in the full four dimensional model, and it is for this reason we feel 
this approach will prove very useful in studying the much more complicated electroweak 
model at high temperature. In this paper we denote the 3d scalars and the 4d scalars 
0. The 4d gauge and scalar couplings are g and A, respectively. 

For Qg^/\ < 0(1), by introducing an auxiliary field x we eliminate the explicit </7^ 
term and develop a systematic 1/A^ expansion for Veff and compute the leading and 
next-to-leading corrections in this expansion. By matching linearly divergent integrals in 
the effective 3d theory with the O(T^) corrections in the 4d theory we show how to obtain 
exactly the correct high temperature result of the full 4d model, at least in the g = case 
which is known [6] . The effective 3d model is sufficient to study the nature of the phase 
transition which occurs at values of (p^ and x much less than T^. In contrast to previous 
studies, wc find no |(/)| term in Veff, nor for Qg'^/X <^ N any significant |0|^ term in Veff. 
In fact, to the order we compute, the model admits only a second order phase transition. 
First order phase transitions arc however found for Qg^/X ~ 0{N) and Qg'^/X ~ 0{N3). 

We mention that the high temperature 3+1 dimensional abelian Higgs model has also 
been investigated in [7] as an effective 3 dimensional theory. However, in that reference, 
the e-expansion method is used to study the nature of the phase transition and it does 
not shed any light on the above questions and problems. 

In section 2 we show how the pure scalar result can be obtained in our approach and 
find agreement with an earlier calculation of this author [6]. In section 3 we extend the 
analysis to the gauged case for 6g'^/ X < 0(1). In section 4 we find the leading contributions 
to Veff for Qg^/X ~ O(A^) and 0(A^3). In section 5 we discuss our results and the phase 
transition and our conclusions appear in section 6. The pure scalar example, although 
previously studied at finite temperature [6] , sets much of the foundation for the analysis 
in section 3, including in particular the vacuum polarization effects. 
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2. (j)^ Model to Subleading Order. 

First consider a model with N real scalars 0* with the tree level 3+1 dimensional 
Lagrangian 

m = ^Wd'ctP - ^(0^ - v'f. (2) 

A is a dimensionless coupling and v is the tree level vev of \4>\. The space-time metric is 
■qi^^ — diag(+, — , — , — ). Furthermore, to avoid problems associated with triviality we will 
assume A > and consider the model as an effective low energy model valid below some 
scale A [10]. 

The systematic 1/N expansion allows us to calculate (1) as a perturbation in 1/A^ near 

= [8,9]. Root [9] has evaluated the leading corrections and given formal expressions 
for the next-to-leading corrections to the zero temperature scalar potential in 4,3,2 and 

1 dimensions. The procedure at finite temperature is very similar. 

By introducing a dimension two auxiliary field % we replace the Lagrangian (2): 
L[4>M = 

= -jidcj)) -v)x- (3) 

The auxiliary field has eliminated the 0^ term; the original form of (3) is easily recovered 
by use of the equation of motion for x- 

To calculate the effective potential K//(0) one proceeds as follows. First, using the 
background field method one computes the effective potential as a function of backgrounds 
of and X- Then, the background of x is eliminated by its equation of motion. 

The systematic expansion first requires expanding the Lagrangian L[0, cr] about 
real backgrounds and x thus [8,9,11]: 

0^^0^ + 4^, x^X+^, (4) 




and then deleting terms linear in the quantum fields 0, x- This procedure defines a 
quantum Lagrangian L[0,x], the sum of whose one-particle irreducible (IPI) diagrams 
give what we call the effective action. Since we are interested in the effective potential we 
assume the backgrounds are space-time constants. We also assume non-negative x- 

At finite temperature the model is formally equivalent to a euclidean field theory 
with one compact dimension. When the backgrounds 0^ and x ^''^^ much below it 
is sufficient to study an effective three dimensional theory, a fact that we exploit below. 
Thus, at sufficiently high temperature, our results for four dimensions should be similar 
to those of a three dimensional euchdean field theory with a dimensionful 0^ couphng [7] . 
Indeed, in three dimensions and zero temperature, the leading 0{N) potential has long 
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been known [8]. It has exactly the same form as the sum of finite temperature superdaisy 
graphs that were computed by Dolan and Jackiw [1] for a four dimensional (p^ theory. An 
important point in our approach is that the effect of introducing an auxiliary field x is to 
shift the mass term and as a result there are no infrared divergences in this formalism. 

To perform the dimensional reduction from 4 to 3 dimensions, one can follow [7] . One 
decomposes the fields as follows, 

4>\xn = Y.^n{S)Mr), X{xn =J2Xn{x)Mr), (5) 

where ipn{T) are a complete set of periodic functions on the circle. For what we are 
interested in only the zero modes n = 0, for which ipoi^) = "00(0), are important in the 
effective 3d model. This is because for n 7^ the fields 0^ pick up nonzero masses of 
O(T^) when the compact dimension is integrated out in the action. Therefore, truncating 
the spectrum to keep only the zero modes and integrating out the compact dimension 
gives the effective 3d euclidean Lagrangian (/3 = 1/T) 

1 rr7 ^ 1 3Ar 2 AT ^ „ 2^ 

-L[0o,Xo] = -^x' + YXi<f> ) 

+ lmA-^ + xWo + 11 + 0>oXo - (6) 

Recall that and x are space-time constants. Defining the three dimensional quantities, 

V>' = ^/P^\ 0' = ^4^l X = X/P, v' = f3v^ (7) 
we finally obtain the 3d quantum Lagrangian 

rr-- M 3A^ 2 ^ / 2 ~2N 

L[V,X] = -^x' + ^xiv'-v') 

+lmA-d' + XW + ^ + ^'^'X 

Here we have dropped the zero mode subscript on xo- 

Before proceeding we would like to stress some important points. As shown in [7], 
if one integrates out the nonzero modes n 7^ at the quantum level rather than just 
truncating the spectrum, there is a finite O(T^) correction to the mass of 0o- This 
correction is very important, in fact it is the term that gives symmetry restoration at 
high enough temperature. We will obtain this term another way, in analogy with what 
one does in effective low energy theories of the strongly interacting standard model or 
effective four dimensional supergravity models inspired by string theory [12]. The three 
dimensional field theory will be divergent. Here we will regulate all the divergent integrals 
by introducing the same scale A. In fact, the only divergent integral that will be important 
is linearly divergent and we will regulate it by simply using a sharp momentum cutoff.|^ We 

"A different scfieme, sucfi as Paufi-ViUars, wifl in general give a different coefficient for tlie 0(A) 
result. However, the precise coefficient will not be important in what follows [12]. 



3A^ 
2A 



(8) 
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then give a physical interpretation to A, i.e. the scale at which the full four dimensional 
physics becomes important. Thus, A is proportional to T. This is in complete analogy 
with, for example, the effective 4d theories where the regulating scale is taken to be 
of order the compactification scale. In our case the identification can be made precise 
because the corresponding 4d correction is well known. We will see that the identification 
A = 7r^T/6 will reproduce exactly the results from the four dimensional integrals. In 
addition, since the four dimensional model (2) is only valid up to a scale A we must require 
A < A, i.e. that T is sufficiently small. Finally we note that the 4d effective potential can 
be obtained from the three dimensional one by dividing by /? and using (7). 

To compute the 3d effective potential we must sum the one-particle irreducible (IPI) 
diagrams of (8). For readabihty, and for later use in the abelian Higgs model, we briefiy 
describe the steps that lead to the next-to-leading result for rather than just stating 
the result. Throughout, Tr stands for momentum and internal space traces. 

Integrating out the 0* in this model gives the 0{N) corrections to Veff to all orders 
in h and defines an effective Lagrangian which only has quantum dependence on x- The 
part of this effective Lagrangian which is quadratic in x includes the tree level part and a 
one-loop vacuum polarization part. The one-loop x corrections give the next-to-leading 
corrections in A^, also to all orders in h. 

Because (8) is quadratic in (f, the 0{(px) mixing term can be removed by making a 
change of field variables [6] . This has the effect of replacing the mixing term by 

+ ^ , _ _ / , (9) 

i-d' + X + x/VN) ^i-d' + x) 

where the last replacement holds to next-to-leading order in the expansion. 

The gaussian integral over (p is now straightforward and its vacuum polarization con- 
tribution to the X propogator at nonzero external momentum is known [8,9] in 3d, 

\ r 1 A 1 _i 

sm 



The final result for the leading + next-to-leading potential is [8,9,6] 



(10) 



Veff = -^X' + ^X{^' - v') + -Trln(-a^ + x) + Vmu (H) 



where the next-to-leading contribution from the gaussian integral over x is 




1 \^ A^V3 



1 - 4x/a2 / + X 



(12) 



Root [9] has shown, without exphcit calculation, that the 3d field theory is renormal- 
izable to this order in A^, as it should be. Although we will give a physical interpretation 
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to our regulating scale, renormalizability of the 3d theory places strong constraints on 
the type of linearly divergent contributions we can obtain. For example, a contribution 
of 0(A<^^) is by itself not renormalizable in our formalism. How this enters will be seen 
in the explicit calculations below. 

Although we did not find a simple expression for Vnu, we found that the dominant 
contribution comes from the large external momentum limit of (10), —d^ ^ x, not the 
zero external momentum limit. Although this seems entirely reasonable in the limit of 
vanishing x, why it should be so for sufficiently large x is not apparent. Indeed, we 
would like to stress the importance of vacuum polarization effects at nonzero external 
momentum. If in (12) we had kept only the zero external momentum part of the % field 
vacuum polarization, 11^^(0) = A/ (487ry^) then we would have obtained an incorrect (and 
nonrenormalizable) answer. 11-^-^(0) when used for the whole momentum integral from 
to A gives a dangerous 0(AA'^/y^) contribution in the expansion of the log in (12). In 
fact, the leading field dependent contributions seem to arise from the most ultraviolet 
divergent field dependent term in the high momentum expansion of H^^. One can check 
this assertion when x is sufficiently large by the following method. In this limit, the 
arcsine term in (12) is small for all values of the momentum, and we expand the In using 
ln[l + x] = X + .... The arcsine possesses two different expansions depending on whether 
—(9^ is bigger or smaller than 4%, and accordingly the momentum integral implied by the 
Tr must be broken up into two different regions. The different contributions can then be 
evaluated with the result that very little error is made is using only the large momentum 
limit of (10) everywhere. This effect will also be very important in the next section. 

In the large momentum limit, we have 



(-9^ + x)nxx(-9^) 

We then write eq. (12) as 



AV-52 A^X 



48 



V, 



Ntl 



-Trln 
2 



A 



A 



127r 



- -Trln 
2 



(13) 



(14) 



m 



It can now be argued that if A/48 -C A and x is sufficiently large the 0{'\J —d"^) term 
the log can be dropped [13]. This is what we will assume. 

We can now give an explicit expression for V^//, up to the approximations we have 
made by using the following result for three dimensions (x > 0): 



Trln(-a^ + X) 



d^k , -2 . Ax 
ln(A;2 + x) 



3 

X=^ 



+ const. 



(15) 



(27r)3 ' 27r2 Gtt 

A sharp momentum cutoff was used to evaluate the field dependent divergent part. Adding 
everything up we obtain, up to a constant. 



3iV 2 A^ / 2 
e// = -^X + yX(^ 



N 



Ax 

47r2 127r 



X^ 



AA 
12^ 



47r , 
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12n , 



3 



• (16) 



To write the effective potential for ip"^ alone we use the equations of motion for x to 
eliminate it. To the order we are working it is consistent to use the leading order equation 
for X to eliminate it [9]. dVejf/dx = from the terms proportional to gives: 



If we rewrite the 0(1) terms in Ve/f using this equation, and ignore A/(127r) in com- 
parison with y/x then we obtain up to a constant, 



eff 



3N ^ . 2 



+ iN + 2) 



Att"^ 127r 



/ 



This is the leading + next-to-leading result for Ve// for a pure scalar theory in 3d 
when y/x is large enough. To obtain the high temperature four dimensional result, we 
recall the high temperature 4d leading order mass-gap equation [1,6]: 

Using (7), we see that the identification 

A = 7r^T/6 (20) 

in (|T^ reproduces exactly the four dimensional result (0). In fact, with this identification, 
(|l^) reproduces exactly the leading and next-to-leading 4d finite temperature scalar 
potential found in [6]. The computations here have been much easier. This potential was 
studied in [6] and admits no first order phase transition for the range of 0^ and T where 
it is valid (i.e. not too close to the origin when the temperature is near T2). Similar 
conclusions appear in [7,14]. Our result is essentially due to the fact that Vefj contains 
no O(|0p) terms for < 3x- 

The potential (18) is valid when A/(127r) is ignorable compared to ^/x- We will give a 
more extensive discussion of the phase transition in the more general gauged case, in the 
last section. 
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3. Abelian Higgs Model, ^g^/X < 0(1). 

Now consider the 3+1 dimensional gauge invariant scalar QED Lagrangian 

L[0, A] = ^6,,d,<l^^d^<l^ " ir^'^' - - l^'^'^F,. - ge,,{d^ct>^WA^ + \g^^A\ (21) 

Here, g is the gauge coupling, i — 1...N, are real as before and e^j is antisymmetric 
with the nonzero components 612, 621, €34, 643, ... all having magnitude 1. The field strength 
and covariant derivative Djj^ are given in terms of the gauge field A'^ by 

F^^u = d^A,-d,A^. (22) 

This model is nothing but N/2 copies of the simple scalar QED Lagrangian discussed in 
many textbooks, and also in ref. [1], but with only one gauge vector field. 

We would like to calculate the leading + next-to-leading corrections to Veff in the 
1/A^ expansion at high temperature. As in section 2, we will do this by first writing an 
effective 3d theory. This three dimensional theory will turn out to be an abelian Higgs 
model with an extra scalar. This can be understood intuitively because in 4 dimensions 
the gauge field has three massive components, while in three dimensions it only has two. 
Hence, to get the same light degrees of freedom we need an additional scalar. 

To proceed, we introduce an auxiliary field, as in eq. (3), to obtain 

L[^,X,A] = ^(d<pf + ^x'-^i<P'-v')x-\F'''F,.-ge,,{d,^^)^ (23) 

We then expand this using (4), and delete all terms linear in and x. Since we are 
interested in the scalar potential we need not keep a background for the gauge fields. 

We also need to gauge fix. We add the gauge fixing term 

Lg.f. = -^{d^A^^ + age,,4>'<P'f. (24) 
Za 

Here, a is an arbitrary parameter. Although the calculations simplify in the Landau 
gauge a — > 0, this extra parameter is useful in checking that physical quantities such 
as the critical temperature are gauge fixing independent. The ghost term for this gauge 
fixing is [1] 

Lgh^9{d'' + agY)9, (25) 
where 9, 9 are Grassmanian ghost fields. We assume non-negative g'^, a and a < 0{1) [4]. 
After performing all these operations, the total quantum Lagrangian is given by 

mx,A,9] = L[lx] + lA,d'A'' + l{l--){dAy 

Z Z a 

+^/(0^ + 2<^>i^y^/N + 4>^IN)A^ - -^eij{d,^^)^A^ 
-^ay2^>^6,,-6«0^0' + 9{d' + ag^(t>')0. (26) 
8 



We dropped all total divergences and assumed the backgrounds are space-time constants. 
There is no 0{(j)A) term because of our gauge choice. Finally, L[0, x] is the pure scalar 
part of the quantum Lagrangian. 

To write the 3d effective field theory for the high temperature limit of this model we 
follow exactly the approach described in section two. We simply truncate the spectrum 
to keep the zero modes of the compact dimension of the euclidean theory and then give a 
physical interpretation to the scale used to regulate divergent integrals. This truncation 
means do —>■ 0. Then, with the definitions (7) as well as (a=l,2,3) 

^'' = V^A^ p = t^A', ^=^pe, ~g = gl{~^. (27) 
we write the effective 3d quatum Lagrangian as 

LYp.x.ApA = L[0,x]~lAad'A''-l-il--){dA)'-l-pd'p 

I la I 

+ ^a^2^V'eueH<^V' + ^(-^' + V')^- (28) 
L\(p, x] is given by eq. (8), and Aa = A'^, da = 3°". 

Eq. (|28| ) is nothing but a three dimensional abelian Higgs model with an extra massless 
scalar field p. It has a simple canonical kinetic term and its tree level coupling to the 
scalars is simply +g'^ip'^p^/2. We now compute the IPI diagrams of this 3d Lagrangian to 
leading and next-to-leading order in the expansion. It is clear that the gauge sector 
contributes only at next-to-leading order; there are no Feynman diagrams involving the 
gauge fields that contribute at 0{N). The correction from the 3d ghosts involves only a 
simple one-loop calculation. For the rest, we can either first integrate out the A and p 
at the quantum level, or integrate out the 0. We have done both, and present only the 
latter computation. 

To integrate out the ip we must first make a field redefinition to eliminate terms in 
(28) which are linear in 0. For ^ = this is just the replacement (9) to the required 
order. For ^ 7^ we have to shift in such a way as to also eliminate the 0{g^) linear 
terms in 0. When this is done we found no contributions from the O{g'^0^) terms that 
were important at next-to-leading order. Another way of saying this is that there are no 
Feynman diagrams involving the O{g'^0'^) terms that contribute at next-to-leading order. 
Hence we just forget these terms in the Lagrangian. 

In the pure scalar case we saw that integrating out the scalars gave the leading order 
result and generated a vacuum polarization term 11^^ for Here, something similar 
happens. In addition to H^^ we will generate a vacuum polarization Upp for p and a 
vacuum polarization matrix H"^ for the 3d gauge fields. This result can be derived very 
simply as follows. The scalar loop integral in (28) contributes a piece 



-Tr In 
2 



7^2 



(29) 
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where the inverse scalar propogator is given by the matrix 



A 



(30) 



One can expand ( pQ]) as a power series in A\ the terms hnear in A in this expansion vanish 
because tre^j = and only the terms quadratic a A and p contribute at 0(1) in the 1/N 
expansion using trej^ej^ 
gives to 0(1) 



— A^. At zero external momentum the expansion is easy and 



^Tr ln(A-.^ + 5,,x/^) + ^^'(^' + p')Tr(-^ + x)"^ + fAaA,TTd'^d'/{-B' + xf- (31) 

The A^, terms are clearly vacuum polarization corrections to the tree level kinetic terms 
of these fields. The vacuum polarization for x is contained in the first term above, and 
is in fact the same as in the pure scalar case at 0(1). Since we have seen that vacuum 
polarization effects at nonzero external momentum are very important even in the pure 
scalar case we will also keep them here. This can be done by carefully accounting for 
the 3-space dependence of all the quantum fields in (|29|) and the result for the full A, p 
vacuum polarizations is nothing other than the one-loop result given for scalar QED in 
many standard textbooks [15]. Altogether, one has 



o 



-O- 



A 



1 



Q W + xEp + ky + x] 
1 



+ x' 



, ^^^^^ I 



6' 



ab 



l{2p + kY{2p + ky 

f + X W + xW+ky + x\l 



(32) 



Veff then is given by the tree-level piece, a contribution from the scalars loops, 
^TrlnA"-^, a simple one-loop contribution from ghosts ^, and the contributions from 
quadratic integrals over A, p, x with propogators modified by the vacuum polarization 
effects (p2D. The result is: 

Veff = -^x' + jX{v'-v') + lTTlnA-^ + VMu 

-Tr ln(-^ + a~gV) + ^Tr ln(-92 + + ^pp) + ^Tr ln(A-,i + n,,).(33) 

The tree-level gauge kinetic term is given by 

A;,^ = -SabB" + (1 - a-')dadb + ~g\^5ab. (34) 

The contribution from x is the same as in the pure scalar case, eq. (12), because 11^^ did 
not change to 0(1) in 1/A^. 
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A few remarks are in order: 

• In the Landau gauge a 0, the 3d gauge contributions are given only by the last term 
in (33). In this gauge, the p and Aa contributions are similar to the 4d gauge loop ring sum 
worked out in refs. [2,4]. As mentioned in the introduction, their procedure overcounts 
the 0{h^) graphs. For the g dependent 0{h^) graphs the overcounting is due to an 
additional contribution from the scalars because for the scalar loop contribution they use 
propogators with g dependent vacuum polarization effects. Our method of computing Veff 
automatically avoids this overcounting. In addition our contribution from (33) includes 
much more than the ring sum in [2,4] because our vacuum polarization effects depend 
on X which is a self consistent solution of dV^ff/dx = to all orders in U, whereas the 
vacuum polarization used in [2,4] is of 0{fi) only. 

• In a general gauge, a 7^ 0, the contributions to \4// due to nonzero Yiab are independent 
of a. This is a result due to gauge invariance, i.e. k"''n.ab{k) = formally holds for (32) - 
as it also should for a properly regulated Hab- We have A";,^ + Hab = Sab{—d^ + g'^^p'^) + 
[(1 — a^^)dadb + Hah]. The log of this can formally be expanded as a power series in the 
term in square brackets. Since Hab is transverse, tr[(l — a~^)dadi, + Hab]" factorizes so 
that there are no mixed a-dependent and 11 dependent terms. Altogether, this means 

Ttln[Aj + Uab] = TrlnA-,i + 2Trln(-92+^V' + n)-2Trln(-a^ + ^V') 

= Trln(-a^ + a^V^) + 2Trln(-^ + ^V^ + n). (35) 

Here, Hab = {^ab — kakh/k'^)H which is true due to ¥'Tlab = and 3d euclidean invariance. 

• Renormalizability of the model places strong contraints on the type of linearly divergent 
corrections to Veff- The reason is the same as in the pure scalar case and as in section 2 
we cannot obtain, by itself, a 0{\Kip'^) contribution to K//. 

• As already described we will give a physical interpretation to the regulating scale for 
divergent integrals. In general, one should introduce arbitrary parameters for the different 
divergent terms [12]. In the pure scalar case we did not do this because it turned out 
to be unnecessary. In our computations for this section it turned out that the match 
between the hnearly divergent 3d integrals and the finite corrections of the full theory 
required that the hnearly divergent part of IIpp be regulated not with A — > 7r^T/6 but 
with A' 7r^T/3. In general, we should also use unspecified sharp momentum cutoffs 
for the linear divergences from the final p and Aa loop contributions. However, the same 
cutoff A — > that gave the O(T^) corrections in the scalar case also worked here. 
This should be kept in mind in what follows. 

Many of the traces in (34) are straightforward to compute. Using (15) we find up to 
field independent terms, 

[X + ag (f )2 -x^\ , 

Trln(-a +ag^) = ^-^ 



TV In A 



N 



Ax 
27r2 



3 



+ 



ag^Aif^ 
27r2 



1 
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Trln(-(9 + g (f + IIpp) = _ „ ^ —. (36) 



3 



'^^^ 27r2 67r 



Notice that the contribution from the 3d scalar p occurs only in the combination ip"^ — 
y^/47r, a fact that is demanded by multiplicative renormalizability of the 3d theory.^ 

We cannot compute the last term in (33) without first evaluating Uab at nonzero 
external momentum. This can be explicitly done by following [15]. Gauge invariance 
forbids a linearly divergent gauge mass correction and in dimensional regularization the 
vacuum polarization displays no poles at ci = 3[|. Since Uab = {^ab — kakb/k^)Il, we simply 
compute n = IS^^Uab- We find 



-Vx + — H~ 



(37) 



P \^Jl + Ax/k\ 

n = 0, (3^ 



This has the zero momentum value 
and the high momentum limit 

n ^ + (39) 

47r 32 ^ ^ 

In the zero momentum limit Uab = 00? and if this is wrongly used in the whole momentum 
integral in (33) we will obtain a nonrenormalizable answer for the 3d potential. This is 
most easily seen in the Landau gauge for which all linearly divergent corrections from 
the 3d scalar sector come in a combination proportional to — The 3d gauge 
loop contributions for Uab = however contribute 0{g'^Aif'^) corrections which are not 
renormalizable. Thus, keeping only the zero momentum limit of Ilafe is incorrect. 

We believe that the dominant contributions to K// from the 3d gauge loop integral 
arise from the high momentum limit of Uab', we justify this as follows. For g'^/12'n' ignorable 
in comparison to it follows that U{—d^) is small compared to —d^ for all values of 
—d^. We then expand Trln(— 9^ + g'^ip'^ + 11) as a power series in 11. The momentum 
integrals implied by the Tr must be broken up into the regions —d"^ > 4% and —d^ < 4% 
in accordance with the different expansions 11 possesses in these two regions. At 0(11) the 
different contributions can be computed with the result that only a small error is made 
in using (39) for the full momentum integral [17]. 

''Note that due to our aim of calculating the high temperature Vef / , we identify p with the longitudinal 
4d gauge field so that a bare mass for p is neither necessary nor allowed. 

'^In 3d this last statement can be understood as the requirement that there should be no nonlocal 
ultraviolet counterterms [16]. 

''This agrees with the results given for the transverse components of the photon in [5] . 



12 



When A S> — .JxI^t^ ^ (7^/(64)^ [13] we find, using eq. (35), a sharp momentum 
cutoff A and ignoring constants. 



Trln(A-,i + n,,) 



+ 



~9' {v' 



71^ 



Sir 



(40) 



Adding everything up, we obtain for the 3d K//, 
3AA „ „ „ Av 

Veff = 



-^X +yX(^ -.) + ^(Ar + 2 + ^, 




3 

3x + ^(5'-A/27r2)y + (Ar-l)xt 



127r 



X + 



Xv 



r,2 



+ 



6 1277' 



:(A'-A) +2 x + 



AC 



r,2 




+Va + neglected terms + 0(1/A^) terms. 
Va is the pure gauge fixing dependent part. 



1 



127r 



(41) 



(42) 



and satisfies V^(x = 0) = Vaif"^ = 0) = 0. The "neglected terms" in (41) refers to the 
(small) errors that were made in only keeping the large external momentum limit of the 
vacuum polarizations, as well as other approximations. Here, we also used the leading 
order result for x, eq. (17), to rewrite the next-to-leading terms in Veff - a substitution 
that only changes the 0(1/A^) terms. 

To get the high temperature 4d Veff from (41) we must match the A, A' terms with 
the 4d result for finite pieces as well as multiplying the 3d potential by T and using 
eqs. (7), (27). For the case ^ = we already found A = reproduces the correct 4d 

result. In fact, this identification also works in the g case. This can be checked by 
examining the 0{h) (i.e. one loop in usual perturbation theory) result for which the third 
term in (41) reduces to [6] Alp^{NX + 2A + 18^^^24712. The identification A = tt^T/G 
reproduces the one-loop result for scalar QED given by Dolan and Jackiw [1]. Finally, 
A' = tt'^T/3 reproduces the 0{T'^) scalar QED longitudinal gauge boson polarization given 
in refs. [2,4,5][|. To arrive at (41) we made the assumptions A ^ ^Jx > XlYlix ,g^ jVlix , 
A » - v^/4vr » W(64)2 [^g]^ 



^To compare we must rescale Lp^ — > <f^ /N , cf^N and set N ^ 2. 
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4. Abelian Higgs Model, Qg^/X - 0{N) and 0{Ni). 

The case 6g'^/X ~ 0{N). In section 3 we considered A fixed as N increases. In this 
case we must assume (XN) fixed as N increases. For a weakly coupled model this means 
the scalar self-coupling is very weak. First let us consider the pure scalar case of section 
2. We again work in 3d. To obtain the leading corrections we make the replacements 
A = X'/N, X X/N, X X/N in (8) and (11). Then Vmi is at most of 0(l/iV); the 
remainder of (11) is of 0(1) and 0(l/\fN). Keeping only the 0(1) and 0{l\fN) terms 
the answer can be exactly computed. For non-negative x, (H) gives: 

V,^f = + Ixiv" - ^5') + ^ ^ + const. (43) 

X should be eliminated by its leading order equation from the the first three terms above. 
For 5^ 7^ the computation of all 0(1/V^) terms is nontrivial. The replacements A = 
X! /N, X ~^ xl^ ill (33) does not give all 0{1/\/N) terms. The 0(1) terms can however 
be found in this way. The final expression is 

- f(^V27r' + <^')i + 2\ip\^\ + const. (44) 

X is given by the same solution as in the ^ = case. Eq. (44) is a independent to 
this order and is valid for non-negative Xi s/X ^"^^ ^ ^ Xj^'^ a^id (p'^ ::$> ^2/(64)^. The 
matching conditions for the linear divergences is the same as in section 3. 

The case 6g'^/X ~ O(A^i). Here we must assume (AA^i) fixed as N increases. We can 
compute the leading corrections by using the replacements A = X'/Ns, x xl^^ ^"^^ 
X — > x/Ni in the results of section 2 and 3. For the pure scalar case, V^ti in eq. (11) gives 
at most 0{N~3) corrections. The rest of eq. (11) gives 0{N^) and 0(1) corrections. For 
g ^ 0, the gauge-loop ring sum included in (33) gives 0(1) and also more subdominant 
corrections. Thus the 0{N^)+0{1) corrections can then be easily found in this way. The 
final expression is 



1 ,^1 ,^1. 3 

2 



3iV3 m 2 mAx X 
^eff = -rrriT^ + ^xiv - v) + 



2mX 2 ' 47r2 127r 

+^^' - ^ f(^V27r' + ^')t + 2\^\^] + const. (45) 

X is found from its 0{N^) equation of motion, and the comments after (44) apply. 
For negative x have to give meaning to eq. (15). One possibility is to simply use 
Trln(— 9^ + x) = Ax/27r^+const. This leads to a real potential, but we have no rigorous 
way of justifying this prescription. 
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5. The Phase Transition. 

We first discuss our results generally before investigating the phase transition. As 
mentioned in the introduction, ordinary 3+ Id perturbation theory in h is not reliable 
for high temperature studies of spontaneously broken theories because it suffers from 
infrared divergence problems just near the value of the temperature at which the phase 
transition from the symmetric phase to the broken symmetry phase occurs. One can try 
to circumvent this problem either by trying to resum certain infinite classes of Feynman 
diagrams to all orders in h or by performing perturbation theory in another parameter, say 
In order to extract reliable results such a new perturbation series, if it is calculable 
to any given order in 1/N, must be renormalizable and avoid infrared divergence problems 
to any given order in 1/N. 

For a weakly coupled 3+ld abelian Higgs model at high T with 6g^/X < 0(1) we 
introduced an auxiliary field x the Lagrangian so as to develop a systematic expansion 
for Veff in 1/N, where N is the number of real scalar fields. This was the content of section 
three. The 0{N) and 0(1) corrections in this expansion are renormalizable and avoid the 
infrared divergence problems of perturbation theory in h. The result can be expanded 
to all orders in h and in fact corresponds to summing certain infinite classes of ordinary 
Feynman diagrams in h. The 0(N) result for Veff incorporates the sum of all superdaisy 
graphs of Dolan and Jackiw [1] at 0(N); the next-to-leading order result incorporates 
much more, i.e. the sum of all Feynman diagrams of 0(1) as well as 0(1) contributions 
from the superdaisies. The 0(1) gauge contributions are given by a gauge loop ring sum 
(but not quite the simple one of [4]). The beauty of our systematic approach over trying to 
add up some infinite classes of diagrams by hand is that there is no overcounting of 0(/i^) 
- or any other - Feynman diagrams and in addition the contraint of renormalizability helps 
to locate all relevant contributions to Veff in a consistent computation. 

For weak A and strong g it is not possible to define an expansion in 1/A^, and for a 
weakly coupled theory it is only possible to set Qg'^/X ~ 0{N) if the scalar self-coupling is 
not just weak but very weak. In this case the expansion of section 3 is unreliable because 
the dominant contribution is no longer from the pure scalar fields; this is demonstrated 
soon. For a very weakly coupled 4d pure scalar field theory at high temperature with 
XN ~ 0(1) it is possible to write Veff as the sum of 0(1) and 0(1 /^/N) terms, with 
subdominant terms of 0(1/A^). In section 4 these were found to give just the daisy sum 
(not superdaisy sum) of Dolan and Jackiw [1]. If computed consistently to 0{1/\/N) the 
daisy sum cannot lead to spurious field dependent terms that come from overcounting of 
Feynman diagrams in h. The daisy sum is of course renormalizable and avoids infrared 
divergence problems. In the gauged case with very weak A and 6g'^/ X ~ 0{N) a calculation 
of Veff to 0(l/\/iV) is complicated because the contribution from the gauge sector involves 
more than just the gauge-loop ring sum of section 3. In fact, if one just adds the pure 
scalar daisy sum result to 0{l/\fN) with the gauge-loop ring sum to 0(1/ VN) then the 
result is multiplicatively nonrenormalizable in the context of the 3d field theory. We did 
not find a simple expression for the 0(1/ VN) corrections involving gauge fields and do not 
know if these corrections have infrared divergence problems or not. The 0(1) corrections 
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however have no such problem and in addition give a gauge fixing independent effective 
potential to leading order. 

For Qg^ /\ ~ 0{Ni) it is possible to easily determine the 0{N^) and 0(1) contribu- 
tions to Ve//, with subdominant contributions of 0{N^^). These were found in section 
4 and have the interpretation of the daisy sum of the pure scalar case plus a gauge-loop 
ring sum taken in the limit of vanishing (field dependent) scalar mass. The effective po- 
tential to 0(1) has no infrared divergence problems, is renormalizable and is gauge fixing 
parameter independent. 

We now study the implications of our results for the phase transition in the weakly 
coupled abelian Higgs model. As noted in the introduction, we derived our results by 
working with an effective 3d field theory (given by cq (28)) and matching linearly divergent 
terms in the effective theory with finite O(T^) corrections in the full 4d model. The 
"neglected terms" therefore incorporate subleading terms that were dropped due to the 
approximations made in calculating the 3d effective potential, as well as finite terms that 
the 3d effective theory cannot account for and that arise when the full 4d model is used. 
We stress that when x <^ and (jP' <^ T^, the 3d effective theory is sufficient to calculate 
the leading (i.e. T dependent) finite corrections. This is of course the region of interest 
for the phase transition. 

The case &g'^/\ < 0(1). Our main result is the leading + next-to-leading order high 
temperature Veff{(j)'^) for a 3-|-ld abehan Higgs model with N real scalars, computed in 
the systematic 1/N expansion. For the tree Lagrangian given by (21) we obtained 



eff 



^ ( (sx + ^{v'- rVl2)] \{N- l)x^ 



TjQgVX) 
127r 



X + ^{v' + rVl2) j \2(x + ^{v'- T'/U) 



+Va+ neglected terms -|- 0(1/A^) terms. 



(46) 



where K, is 



T 



127r 



ix + ag )2 -x^ -a^g 



(47) 



and x(0) is the solution of the equation dVeff/dx — 0. To next-to-leading order it is 
sufficient [9] to use not the full solution but the 0{N) solution of (19), 



XT 
487r 



1 + 



327r2 /1202 I2v^ 



X 



2^2 



J^2 



+ 1-1 



(48) 



There is another solution of y/x which comes with an overall minus sign, but this is 
unphysical [6,8,9]. Our Veff is valid for T > ^ > XT/Un, g^T/12n and ^ 0^ _ 
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T^/Att > 5(2T2/(64)2. Also, eq. (48) was used to rewrite the 0(1) terms in (46), 
which adds \{dVi/ dxf I'iN to the 0{1/N) terms, where Vi is the 0(1) part of (46). In 
the Landau gauge, and near T = \/V2v, this particular 0{1/N) term is negligible for 
X/T2 > (A + V)(l + (6(?VA)i/2)/37V2_ 

In our approach vacuum polarization effects at large external momentum played an 
important role in the same way as they did for the previously studied [6] = case. 

It is known [6] that at order iV, the potential (46) admits no 1st order phase transition. 
The 0{N) potential is exact in the small x hmit (up to 4d corrections). To 0{N), 
dV^f f/d(f)^ = dVeff/d(f)^ = Nx/2. At 0^ = this vanishes at T| = 12v^. For T > Vl2v 
the origin is a global minimum, and at T = T2 the potential grows away from the origin. 
For T > T2 the point % = is away from the origin and this has the interpretation [8,9] 
as the symmetry breaking minimum below T2. For the case g = 0, Root [9] has shown 
for the 3d case that the point x — ^ remains a local minimum at next-to-leading order. 
This analysis did not require a computation of 14//. Root examined dV/dcfP' in the limit 
of vanishing x ^"^^ showed that the leading order gap-equation for x iy^ our case eqs. 
(17) and (19)) was sufficient to show that \4///<^0^ = still has a solution at x = at 
next-to-leading order. His analysis can be extended to our gauged case. We will not 
present a detailed analysis here but instead refer the reader to Root and also ref. [18] 
where, following Root, an analysis of the small x limit in the full 4d high temperature 
abelian Higgs model has been performed with the result that X — remains a point of 
vanishing dVeff/dcj)^ at next-to-leading order. 

We believe the fact that x = remains a point of vanishing dVeff/dcj)^ is however 
insufficient to deduce a second order phase transition to this order. As mentioned by 
Root, for sufficiently large N the 1/N corrections cannot overwhelm the leading order 
conclusion of a second order phase transition. For N not arbitrarily large we would like 
to know the global properties oiV^ff away from the point x = 0; ^^'^ iii particular if there 
is a point away from the origin at T = vT2f where the 1/A^ corrections can produce 
a new minimum and possibly even result in the breakdown of the 1/A^ expansion. In 
addition, for T > \/V2v the point x = never occurs, so the analysis of [9] is by itself 
insufficient to deduce 0^ = remains a global minimum for N not arbitrarily large. Our 
computation of Ve// sufficiently away from x = gives global information that Root's 
analysis cannot give. In addition, (46) appears to characterize the correct behavior in the 
limit of vanishing x- Assuming no pathological behaviour occurs in an exact computation 
of Ve// at next-to-leading order for ^Jx < \T / 127i , g'^T / 12n we might expect that (46) 
gives a good description all the way down to ^/x — ^^'^ ^Iso 0^ = at ^/x — 0. 

With these points in mind we now investigate if, more generally, (46) can exhibit a 
first order phase transition for 6g'^/X ~ 0(1). 

The critical temperature T2 is given by the vanishing of dV/dcj)'^ at 0^ = 0. Let us 
write Veff — Vn + Vi + V^, where Vn is the 0{N) potential. Then since dV^/dx — we 
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obtain 
dV 




0(1)'' 

T y 



+ 



Nx 



T 
Stt 



(49) 



As indicated by our preliminary remarks, this vanishes at x = and at the origin this 
translates to the gauge fixing independent result 



(50) 



This is in fact the leading order result again (i.e. as if Vi were absent). The reason is that 
in a general model the 0{1/N) terms in are important in determining the 0{1/N) 
corrections to T2 (this is related to the fact that x = remains a point of vanishing 
dVeff/dcj)'' at 0(1)). 

For T > T2, 0^ = is the local minimum. At 0^ = 0, V„ = so K// is gauge 
parameter independent at this minimum. For there to be a first order phase transition, 
eq. (49) must have zeros for T > T2 and, at T = T2, dV/dcp'^ should be negative away 
from 0^ = 0. Let us work in the Landau gauge, a — >^ (in an arbitrary gauge a should 
be not more than 0(1) [4]). Then at T = Ts, 



dx 12 I ^ ' 



An 



(6^VA)) - g VX (3^3 - 1 + 2(6^VA)i) . (51) 



X increases with 0^, hence a necessary condition for dVeff/d(j)'' to become negative at 
T — T2 is that (51) becomes negative. At (7 = 0, (51) vanishes only when ^/x ~ nT2/6 
by which time the effective theory is no longer valid and the positive contribution to 
dV/d(p^ from the 0{N) term in (49) already dominates. Hence for g — no first order 
phase transition is possible by T = T2. For Qg'^/X ~ 0(1), g' <^ 2n, the same conclusion 
can be reached, namely that within a consistent 1/N approximation the leading+next- 
to- leading Vgff for our weakly coupled abelian Higgs model does not admit a first order 
phase transition. (Note the absence of an O(|0p) term at T = T2 - see eq. (48).) 

If Qg'^/X is of 0{N) then (51) can become strongly negative, indicating the possible 
breakdown of the 1/N expansion. Note also that the "ncxt-to-leading" corrections to 
the coefficient of x in (46) actually compete with the 0{N) coefficient and therefore our 
results cannot be used in this case. 

The case Qg'^/X ~ 0{N). Our main result is the high T contribution to \4// which 
survives in the limit of arbitrarily large N. For the tree Lagrangian given by (21) we 
obtained (after eliminating x): 



eff 



:^0^(0^-2^^ + r76) + ^0^ 



g^ 
'l2n 



{Ty6 + r)2 + 2\(j)\-^ + O(lViV) + neglected terms. 



(52) 
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Note that for weak coupling one must keep g and (AA^) fixed as increases. As mentioned, 
the potential above is a independent to 0(1) in the way we computed it. One can also 
check that it is gauge fixing parameter independent with the slightly different gauge 
fixing considered in [4]. Eq. ( ^2]) is strictly valid for 3> 0^ ^ T^(yf^/(64)^. The critical 
temperature T2 determined from (|5^) is given by 



12v' 



1 + 3(6//AAr) 



rV2 



(6//AAr) 



(53) 



We assume (7 <^ 27r so that T2 here is lower than the leading order critical temperature 
for the 6g^/\ < 0{N) case, vT2f.Q At T = T2 the minimum is no longer at the origin 
and occurs at the point 

^ = IgiQgVm, (54) 

indicating a first order phase transition whose strength grows with the gauge coupling. 

Since we did not compute subleading corrections to K// all we can say is that the 
model exhibits a first order phase transition for sufficiently large A^. This is in contrast 
to our results in the 6g'^/XN < 0(1) case. 

The case Qg^X ~ 0(A^3). This case is the most interesting because, unlike the 
case 6g^/X ~ 0{N), the next-to-leading corrections are easily determined. The high 
temperature potential is 



4! ^ 
_^ 
127r 



3 

127r 



oix) 



+ 0^)^ + 2|0|^ + 0(A^"3) + neglected terms. 



Here, 



X 



xnI 



(55) 



(56) 



We note that for weak coupling, AA^s must be held fixed as A^ increases and that the result 
for Veff is gauge fixing parameter independent to 0(1). Finally, note the appearence of 
the step function 9{x) = 1 for x > , otherwise, in the expression for K//- As 
described in section 4 there is an uncertainty involved in defining the gaussian integral 
over the quantum scalar fields when x is negative; we have chosen an answer that in light 
of the results for the cases Qg^/X < 0(1) and Qg^/X ~ 0{N) gives a physically appealing 

•^At the origin and at T = T2 the auxiUary field x that was efiminated from (44) is negative, and one 
may wonder if (44) and (52) are well defined in such a situation. However, the leading divergent part 
of Ti{—d^ + x)~^ is independent of x and well defined. Hence, the linearly divergent term in (44) is 
unambiguous for negative x- In fact, to this order the quantum contribution from only the scalars is just 
the 0{T^) shift to the tree-level mass of the scalars already derived by a one-loop 0{h) calculation in 
[!]• 
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interpretation. Therefore, our results here are somewhat speculative. We assume g <^2ti 
in what follows. Also V^ff is strictly valid for > ^ ^nd (jP' > g'^T'^ / {Q^Y . 

The leading order critical temperature is the same as in the case Qg^/X < 0(1), i.e. 
T2 = \/l2v. The 0(A^3), potential given by the first term in (55), does not admit a first 
order phase transition. The next-to-leading order T2 is modified. We find the solution 
dVeff/dcj)'^ at 0^ = occurs at negative x and it is again given by relation (|53| ) (but note 
that it reduces to 12f^ = T2 for N 00 here). At T = T2 the global minimum is away 
from the origin and at (^). For Qg^/X ~ 0(A^3) we have 

\(P\min n T\T /cr'yN 

— — > as A* ^ 00, (57) 

indicating a first order transition that gets weaker and weaker as A^ increases. 

We end this section with the following discussion. If we chose a prescription without 
the step function in (55) the potential is complex at or just above T2 given by (53) and 
close enough to the origin. Since we have all terms to 0(1) the imaginary part cannot be 
removed at sufficiently large A^ by subdominant corrections. We do not know the meaning 
of this result; neither do we know in cases a) and b) if subdominant corrections we have 
not computed can lead to the possibility of a complex potential at T just above when the 
phase transition is supposed to take place. However, one can check that in the case here 
even without the step function in (55) the expression for T2 is real to 0{N~s) and given 
by (53) and also that the value of Veff at the minimum (54) is real to 0(1). Our physical 
results, to the order we are working in, are therefore not really so speculative. 



6. Conclusions. 

We have found that for large A^ the weakly coupled abelian Higgs model can exhibit 
either first or second order phase transitions with first order behaviour possible for 6g^/ X > 
0{N). For Qg^/X ~ 0(A^3) the model exhibits a weak first order phase transition at large 
A^. Intuitively, this only says that the gauge coupling must be large enough so that the 
gauge contributions can dominate or at least compete with the pure scalar contributions 
to Veff. For the three cases a) Gg^X < 0(1), b) Qg^/X ~ 0{N) and c) 6g'^/X ~ 0{nI) 
we estimate that our results are reliable whenever the parameter a) A^~^, b) A^~2 and 
c) A^~3 is sufficiently small. An inspection of our results for case a) shows that even for 
A^ = 4 the next to leading corrections to the coefficient of 0^ in Veff compete with the 
leading correction to the coefficient of 0^ in Veff. We do not believe any of our results 
can be directly applied to the A^ = 4 case; however they can give some indication of what 
happens at A^ = 4 in a very weakly coupled theory with (7 -C 1, A <C 1. In this limit, 
cases b) and c) admit only a very weak first order phase transition. For case a) it can be 
argued on the basis of the 3d field theory (and dimensionality of the 3d couplings) that 
the coefficient of 0^ in Veff to lowest order in the couplings and to all orders in 1/A^ is 
given by the one-loop 0{h) graphs. This yields a critical temperature given by [1] 

T|(l + 2/A^ + 18//AA^) = 12v'^. (58) 
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For = 4 and Qg'^jX ~ 1 this gives a critical temperature much below \/Vlv and 
demonstrates that our analysis of case a) is insufficient to rule out a significant first order 
phase transition. Assuming that 0{1/N) corrections to (46) are negligible it is possible 
to use (46) to study the phase transition to temperatures below T = \/Vlv. The leading 
order solution ^Jx of (48) becomes negative and eventually complex near the origin as 
we lower the temperature below T = \fVlv. To avoid this problem one can use the full 
solution for x from (46), not just the leading order one. This approach is not entirely 
consistent but the result of the analysis can be shown to be again a second order phase 
transition to T > T2 given by (|55|). This only demonstrates that any first order behaviour 
is contained in 0(1/A^) corrections to Vejf which we can expect to be small even at = 4. 
Hence our result for all cases is that for a very weakly coupled model the phase transition 
of the exact model, even at = 4, will be for all practical purposes second order. 

For (7 ~ 1 and = 4, we believe that a reliable determination of the nature of the 
phase transition cannot be given by the 0{N) +0(1) corrections in case a), nor by the 
daisy and ring sum results of cases b) and c). In all cases more subleading contributions 
must be computed than those that we have computed or currently exist in the literature. 
These extra contributions may be negligible even for (7 ~ 1 and A^ = 4 but this is not a 
priori clear. However if the exact g = model admits only a second order phase transition 
then we can expect that for g and Qg^/X < 0(1) the exact model will admit at best 
quite weak first order behaviour. 

Our physical results are in all cases gauge fixing parameter independent, they do not 
overcount any Feynman diagrams and do not suffer from any infrared problems. We found 
no terms of 0{T^\(f)\) in Ve// in any of the cases. One can argue on the basis of the 3d 
field theory, dimensionality of the 3d couplings and renormalizability that any such terms 
must be very small in a very weakly coupled model. More generally, any such terms must 
be vanishingly small in the limit of large A^. 
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Note Added. As we completed this work we recieved a paper by J. March-Russell, 
LBL-32540, PUPT-92-1328, which also considers 2nd order phase transitions in gauge 
theories, but using e-expansion techniques. 
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